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We estimate the nuclear fusion rates for the transition 17 1 ff — > 18 F* from the ground and 
excited molecular states and also consider effect of screening of electric charges of nuclei by bound 
electrons in the case of colliding atoms 17 O and 1 H. We show that the fusion rate for the molecular 
state is vanishingly small while the electron screening drastically changes the rate for the collid- 
ing atoms at kinetic energies up to hundreds of keV. The latter may be of consequence for some 
00 ' astrophysical applications. 
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£N1 , It is a matter of common knowledge that the nuclear fusion rate from a molecular state is strongly suppressed by 
Coulomb potential barrier and such fusion events are undetectable. Our interest in the issue is provoked by potential 
possibility of enhancement of nuclear fusion rates from molecular states due to degeneration of energies of initial and 
final states 0. In this paper, it was shown that the transition amplitude for the d 6 Li — ► 8 Be reaction increases 

' due to degeneration of energies, E^hi — 22.2 MeV and Es Be = 22.279 MeV ("nuclear resonance"). As a simple 

, assumption, the following resonance wave function of d 6 Li was chosen: 
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with rj = ZiZ 2 ay / fj,c 2 /2E. The molecular wave function was constructed as a product of the regular Coulomb solution 
Fo(k,r) and an exponentially decreasing function associated with size of the molecule, 
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Normalization factors N rcs and N mo i were not calculated. It was found out that amplitude transition (taken as 
overlap integral of 'J'rcs and ^/ mo i) increases as 7y 3 / 2 e n - 1147ri ) ; instead of the usual Coulomb barrier factor e^ 17 . It was 
proposed that such an enhancement may even be detectable if the normalization factors N Tes and N mo \ are not too 
large. Apparently, this speculation may be generalized to other relevant molecular states. So, we decided to elucidate 
what is actually going on, using a more rigorous approach to the problem. However we do not consider nuclear forces 
explicitely and associated nuclear resonance. We restrict our calculations by electromagnetic forces emphasizing value 
of wave function discribing relative motion of two nuclei in molecule at nuclear scale distances. Effect of nuclear 
resonance leading to the strong enhancement of the fusion rate could be included later in astrophysical S'-factor. To 
be specific, we consider the hydroxyl molecule that is the system of two nuclei, 17 O and 1 H, with nine electrons 
around. This choice is made taking into account the degeneration between the total energy of two nuclei 17 O and 
1 H (5.607 MeV) and the energy of the "united" excited nucleus 18 F* (5.605 MeV). Solving the relevant Hartree- 
Fock-Roothaan (HFR) equations for the system within the method of molecular orbitals and taking into account the 
actual distribution of electrons, we calculate the effective potential between the two nuclei. With this potential, the 
wave functions for bound and scattering states of the nuclei are obtained. Using these wave functions, we estimate 
the probability for the nuclear fusion 17 1 _ff — ► isp* DO t n from the molecular and scattering states. The latter 
case may be of potential interest in studies of thermonuclear reaction in stars. For stars like the Sun, the common 
accepted Debye-Huckel approximation to the calculation of screening effect is, strictly speaking irrelevant, because of 
huge difference between the scattering scale (nuclear size) and Debye radius, Rd (atomic size). This fact is usually 
of no importance considering that the screening due to the plasma electrons is by itself rather small (see Ref. || 
for more details). On the other hand, the low-lying bound electrons (if exist) do screen electric charges of nuclei 
much effectively. Evidently the concentration of electro-neutral atoms inside a star is very small. Nevertheless, some 
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appreciable fraction of ions with low-lying bound electrons (which mainly define the value of the screening) must be 
present in stellar plasma and therefore can affect the common rate of thermonuclear fusions in the stars. Obtained 
results for the specific case of the 17 1 H — > 18 F* reaction indicate very important role of such screening for value 
of fusion rate and we believe that this effect will be common for others nuclear transitions with heavy enough ions. 



II. CALCULATION OF THE EFFECTIVE POTENTIAL BETWEEN TWO NUCLEI 



Adiabatic approach allows us to split up Schrodinger equation for the system into two equations both for electron 
and nuclear wave functions. The approach valids up to kinetic energy of nuclei about few MeV. Total wave function 
presents as (look for instance ||): 

*(ri, . . . , r 9 ; Ri, R 2 ) = u(ri, . . . , r 9 ) • v(R), 
and the equation for wave function of electrons looks: 
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u(ri, . . . , r 9 ) = E(R) ■ u(r ly . . . , r 9 ), 



(1) 



where E{R) depends on the R = |Ri — R 2 |. Wave function describing relative motion of nuclei 17 and 1 H can be 
found from the equation: 
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v(R) = e-v(R), 
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Here r ai - the distance between i - electron and a - nucleus, - the distance between i and j - electrons. Z a and 
M a - are electric charge and mass of a nucleus correspondingly measured in atomic units. Equation ([j]) gives us so 
called molecular therm E(R) serving as effective potential energy for relative motion of the nuclei. With effective 
potential E(R) in hand we can obtain wave function for nuclei solving (||). At first, let's treat with E(R). Wave 
function u(ri, . . . , rg) presents as Slater determinant whose elements are one-electron functions called molecular spin- 
orbitals. Energy of the system reaches its minimum for those of them which satisfy to Hartree-Fock integro-differential 
equations which however hard to compute for molecules. Roothaan suggested to express them as a linear combination 
of known functions Ufe(r) with unknown coefficients (see Ref J3|). These coefficients (formative matrix C) in its turn 
satisfy to nonlinear algebraic system of equations which solves iterativly: 



FC = eSC 



(3) 



Here F is a matrix containing matrix elements of Fock operator F dependent on internuclear distance R. (see appendix 
A). It depends also on unknown yet C. Overlap matrix S =|| (ui\v,k) \\- Diagonal matrix e contains orbital energies as 
the solution of equation (||). Configuration of molecular spin-orbitals describing ground state of the molecule 17 1 iJ 
is (la) 2 (2a) 2 (3a) 2 (Tr + ) 2 ir~ Q. So we choose molecular spin-orbitals as a linear composition of atomic orbitals with 
unknown yet coefficients Cmf- 



ma = C rnl {ls) + C m2 (2S) + C m3 (2p Q ) + C m4 (lh), 
n+ = (2p+), 

TT" = (2 P _), 



(4) 



where m runs from 1 to 4. Atomic orbitals presented here are the solution of Hartree-Fock (HF) equations for oxygen 
atom H and Is hydrogen function: 



Is 
25 

2po 
2p± 

lh 



(P ls /r)Y 0Q (6,^), 
(P2s/r)Y 00 (6,ip), 
(P 2p /r)Y lo (0,¥>), 
(P2 P /r)Y 1± (6,<p), 
e-r'Yooie'^'). 



2 



The coordinate frame centered on the oxigen atom is used. The primed variables refer to the hydrogen nucleus. 
Radial parts of these wave functions are presented in Fig. |l|. All two-centered integrals were evaluated using following 
presentation for hydrogen function B : 

e - r > = j^e^v Pn(cose) (5) 

n=0 



Here r' — \Jr 2 + R? — 2rR cos 6, coefficients a n (r, R) - are linear combination of half- integer order Bessel functions 
(see appendix B). P n (cos6) - are Legendre polynomials. It allow us to integrate over all angles in appearing integrals 
and to reduce them to one- and two dimensioned. Calculation of the kinetic integrals requires also derivation of 
a n (r, R) over r, all these functions were FORTRAN coded and multiply tested. Used expansion (|^) converges rapidly 
with increasing n. Equation (^) has a 6 x 6 block-diagonal structure for spin-orbitals with different symmetries. 
Thus a- orbitals are the solution of 4 x 4 secular equation (Q) and tt- orbitals can be deduced from the normalization 
conditions. Coinciding symmetries of the ground state of hydroxil molecule (7r-symmetry) and that of the united atom 
Ftor (p-symmetry) for R — > allow us to use the configuration of spin-orbitals (^) for all internuclear distances from 
R = to R = oo. For a lot of molecules these symmetries are different and one has to separate region of internuclear 
distances for regions with different spin-orbital configurations and then to concatenate solutions from different regions. 
So in our case we can easely check out if HFR method works well for O H molecule simlpy running R in the equation 
(|). Molecular energy subtracting nuclear repulsion energy was obtained for R € (0,4) a.u. with step 0.1 a.u. and 
spline interpolated. It is presented in Fig. ^|. Calculated molecular energy smoothly approaches to that of united atom 
(Ftor) (Ep ~ ~199Ry) and gives good approximation for energy of isolated atoms: Eo+h — — 149-Ry. Molecular 
energy reaches the minimum for R ~ 2.2 a.u., and ratio of calculated energy E m i n = —149.86 Ry to the experimental 
value —151.56 Ry Q| at this point is 0.988. So in three points where there are experimental results for the potential 
curve we have a good agreement with experiment. 



III. WAVE FUNCTIONS OF THE NUCLEI 



A. Bound state 



Effective potential between the two nuclei shifted for the energy of isolated atoms and measured in Ry is presented 
in Fig. |^. Numerical integration of the radial equation of relative motion of the two nuclei with orbital moment L = 
was accomplished by means of Numerov scheme. Discrete vibrational spectrum of nuclei was found using method of 
"shooting". Table | presents the results. Energy counts from the bottom of the potential curve. Vibrational energy 
spectrum with AE = 0.04 Ry = 0.54 eV is presented here. Value of ^(0) is measured in atomic units. We also (as in 
p|) note significant enhancement of the rates of nuclear fusion for molecular vibrational excited states. Physical reason 
for such behavior one could simply understand looking on Fig. [§ There are shown ground, first excited and 16-th 
excited wave functions. It is seen that the last one is much more near to the R = point. Nonetheless vanishingly 
small value of^>(fS) forbids the possibility for any cold nuclear fusion from the molecular state 17 1 H even for the 
system with degenerate energies. Estimation for the fusion rate gives 
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A = 1.175 x 10 9 ^ x |*(0)| 2 • sec- 1 , 

MeV-barn 1 wl 

where S(E) -is astrophysical S factor. Numerical factor in this formula is given for value of \l/(0) measured in atomic 
units. 



B. Scattering state 



Practically nuclear fusions from scattering state of nuclei are much more interesting. Thermonuclear reactions in 
stars and nuclear fusions in laboratories on the Earth are such examples. Whereat cross sections of these events 
depend on wherever existing electrons which screen electric charges of nuclei. Thereby correct accounting of electron 
screening of electric charges of nuclei is required. It is a simple matter to show that common accepted Debye-Huckel 
approximation to the calculation of screening effect is irrelevant for not very dense stars like the Sun. This approach 
requires electroneutral plasma, i.e. there are many charged particles within a Debye sphere i?£>. However most of the 
contribution to the nuclear reaction rates comes from the internuclear distances shorter then Bohr radius of colliding 
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nuclei. In the core of the Sun scale of Rd is atomic one and scale of the the Bohr radius is nuclear one. Since 
there are no electrons inside a sphere with nuclear size in solar plasma we come to the conclusion that Debye-Huckel 
approximations becomes irrelevant for above considered situations. Recently this question was discussed by Dar and 
Shaviv M in relation with calculation of neutrino fluxes from the Sun. What actually could screen effectively electric 
charge of colliding nuclei are atomic electrons. To learn effect of atomic electron screening we use effective potential 
between two nuclei obtained in this work. Assuming asymptotic wave function describing relative motion with the 
momenta p of the two nuclei as: 

^ *m(pr + g) 
r 

we integrate equation (|J) and compare obtained ^(O)! 2 with |\l/coui(0)| 2 that has a meaning of the scattering wave 
function of the bare nuclei. In our calculations R m i n = 10 a.u. represents nuclear scale. 

|*co ul (0)| 2 = e- 2 "'2^p 2 . (7) 

Here 77 = m* Z\Zij 'aop, ao - is the Bohr hydrogen radius. Fig. [5] demonstrates the solution of numerical integration 
of equation (|J). The ratio ^jf^y^ is shown in Fig. [| For kinetic energy of colliding nuclei (energies are in the 
center-of-mass system) below 1 keV, ^(O) is many orders of magnitude larger then 'I'coui(O), and for medium kinetic 
energy in the Sun (E ~ 1.4 keV) the ratio q7~~7oj — 10 8 - Since Coulomb barrier (~ 1 MeV), most of the contribution 
to the fusion rates comes from collisions with energies larger then medium energy in the core of the Sun. For example, 
for E ~ 10 keV, the ratio is about 9 and it tends to 1 when kinetic energy of colliding atoms reaches energy about 
1 MeV. Thereby atomic electron screening drastically changes nuclear fusion rates and must be taken into account. 
Physical reason of such effect is "electron acceleration". It means that the difference between energy of united atom 
(ftor) and energy of separated atoms (oxygen and hydrogen) adds to the kinetic energy of nuclei accelerating them (an 
extra attraction discribes by an energy increment U). Thus, the nuclei may be considered as tunneling the Coulomb 
barrier at an effective incident energy E e ff = E + U. Usually, the resulting enhancement of the cross section expresses 
as: 

f(E) = * expW^I}, (8) 

at those energies when U <C E. For considered case of 17 1 i? molecule, U ~ 648eU. And electron screening due to 
the exponential dependence (0) at low energies leads to the very strong enhancement of nuclear fusion rates. Good 
qualitative agreement is observed with values computed from the formula (p|) and from our calculations. 



IV. CONCLUSIONS 



• Althogh the excited state rates of molecule 17 1 H are many order of magnitude larger then those for the ground 
states, all are much too small to be detected even for the system with degenerate energies. • Atomic electrons do 
screen effectivelly electric charge of colliding nuclei, drastically changing their fusion rates for low energies. Although 
obtained results are for the specific case of the 17 O x iJ — > i8p* f us i on we believe that this effect will be common for 
others nuclear transitions with heavy enough ions. 
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APPENDIX A: 



• Matrix element of Fock operator F used in this work looks: 
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F km = H km + 2J Di n (2[km\ln] - [kn\lm]), 

l,n=l 
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where 



Din = ^2 Cu*C in , H km =< u k \h\u m >= J dru k *(r)h(r)u m (r), 

i— occup.orb. 

ri i/i [ [ j j u k *{vi)u m (v l )u*{v2)u n {r 2 ) f- , 1. ^ Z a 
[km\ln] =11 dr 1 dr 2 — — p-^ — r—^ — — , h(r l ) = --A* - 2^ — 



ri -r 2 



• Molecular energy finds as: 



E — ^ E 1 + ^ DkmHmk- 
i=occup.orb. k.m 



e % are the diagonal elements of matrix e. 



APPENDIX B: 



• All two-centered integrals were evaluated using expression (|J). 
Denote: 



2 Mr 

/™0) = \l £-I„+Uz), 9n{z) = - J —K n+ i(z), 

TV y ZiZ 



where I n+ i(z) and K n+ i(z) - are modified Bessel functions of half-integer order. 
Let r' — \Jr 2 + R 2 — 2rR cos 9, then (see for instance 0): 



e -Xr' 



— = K 2n + 1 )fn{zi)9n(z 2 )} P„(C0S( 



Ar „=0 

where Zi = Ar<, z 2 = Ar>. And r<, r>- are minimal and maximal distance from r and R correspondingly. 
One can obtain: 

_ OLn(zi,z 2 ) _ ^ + l)f n (z 1 )g n (z 2 ) + zif ni (zi)g n (z 2 ) - z 2 g n+ i(z 2 ) f n (zi) . 
r 

We need also derivation of these functions over r for evaluating kinetic integrals: 
for r < R: 



da. n (z\,z 2 ) 
dr 

for r > R: 

da n (zi,z 2 ) 



dr 



(2n + l)(n + l)f n (zi)g n (z 2 ) + (3n + 4)zi/„ + i(zi).g„(z 2 ) - 

(n + l)-2f2fl , n+l( 2f 2)/n(2i) - Ziz 2 g n+1 (z 2 ) f n+1 (zx) + zlg n (z 2 )f n+2 (zi). 

(2n + l)(n + l)/„(^2).9n(^i) - (3n + 4)« 1 / n (z 2 )5r n+1 (2;i) + 

(n + l)z 2 g„(zi)/„+l(z2) - Zl^SVi+lO^O/n+l^) + 2ig„(zi)/„ +2 (z 2 ). 
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TABLE I. Discrete vibrational energy spectrum and $(0) 



level,N 


energy, Ry 


*(0),a.u. 
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FIG. 1. Radial parts of orbital wave functions of oxygen 
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FIG. 2. Molecular energy subtracting nuclear repulsion energy 
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FIG. 3. Effective potential between two nuclei 
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FIG. 4. Radial parts of the ground, first excited and 16-th excited wave functions 
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FIG. 5. Scattering wave function at R — of colliding nuclei with electron screening (solid line) and for bare nuclei (dashed 
line) 
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FIG. 6. Ratio ' I ' < ^'"' ( - U ' > for different kinetic energy of the nuclei 
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